Lattice Thermodynamics for Ultra-Cold Atoms 
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We measure the temperature of ultra-cold ^'^Rb gases transferred into an optical lattice and 
compare to non-interacting thermodynamics for a combined lattice-parabolic potential. Absolute 
temperature is determined at low temperature by fitting quasimomentum distributions obtained us- 
ing bandmapping, i.e., turning off the lattice potential slowly compared with the bandgap. We show 
that distributions obtained at high temperature employing this technique are not quasimomentum 
distributions through numerical simulations. To overcome this limitation, we extract temperature 
using the in-trap size of the gas. 



Ultra-cold atoms confined in optical lattices are a 
promising system for studying models of strongly corre- 
lated systems relevant to condensed matter physics. Ex- 
periments have been able to observe a quantum phase 
transition from a superfluid to Mott-insulating state for 
bosons H, Q , the superfluid-to-insulator transition for 
bosons [J], super-exchange the cross-over between 
quantum tunneling and thermal activation of phase slips 
Q, reversible depletion of condensate fraction induced 
by fine-grained disorder Q, and evidence for a Mott in- 
sulator of fermions Q . Straightforward interpretation 
of these results in certain cases has been complicated by 
difficulties in measuring temperature related to strong 
interactions and the lattice potential [Toj . 

Temperature has been measured for both bosonic [ll| 
and fermionic atom gases confined in harmonic traps 
by fitting the momentum distribution, obtained after bal- 
listic expansion, to analytic expressions obtained using 
the semi-classical approximation [l^, [ij] ■ Unfortunately, 
this approximation fails in an optical lattice because of 
the rapid spatial variations in the potential [3, [Hi ■ De- 
spite this limitation, progress has been made on mea- 
suring temperature in an optical lattice using alterna- 
tive techniques. One method involves site occupancy 
statistics [ig], which can be straightforwardly related to 
temperature in the atomic (i.e., tunneling energy much 
less than interaction energy) limit [l3|- Another tech- 
nique is to determine temperature using the visibility of 
the momentum distribution [l8| . This method, however, 
has generated some controversy [13, [H, [13, [HI , cannot 
be used at high temperature, loses sensitivity at very 
low temperature, and does not work for fermions. Fi- 
nally, advances have been made employing quasimomen- 
tum distributions to measure temperature. By fitting 
to quasimomentum distributions determined from a re- 
stricted region of a momentum profile, the condensate 
fraction and temperature have been measured for a lat- 
tice in the limit described by the 2D Bose-Hubbard model 
[^. These experiments, however, were carried out at a 
fixed temperature below the critical temperature Tc for 
Bose condensation. In this work, we extend this tech- 



nique over a wide range of temperatures to an ultra-cold 
gas of non-condensed *^Rb atoms confined in a 3D optical 
lattice. 

In contrast to |2] , we obtain quasimomentum distribu- 
tions [H [H, [H through bandmapping [H H^, a pro- 
cedure in which the lattice potential is turned off slowly 
compared with the bandgap before time-of-flight imag- 
ing. We fit these distributions to a semi-classical anal- 
ysis to determine temperature. To explore a range of 
temperatures in the lattice we vary the temperature of 
the atoms before transfer into the lattice, when they are 
confined in a purely harmonic potential. The temper- 
ature determined using this method is compared with 
thermodynamic predictions for non-interacting particles, 
assuming adiabatic transfer from the parabolic potential. 
We show that this procedure is successful at low tempera- 
tures, but fails at sufficiently high temperature such that 
states with high quasimomentum and localized states are 
occupied. Through numerical simulation, we show that 
this failure is related to a breakdown of bandmapping in 
this regime. At high temperature we demonstrate that 
temperature can alternatively be extracted using the in- 
trap density distribution. 

The paper is organized as follows: Section |T] reviews 
solutions to the Schrodinger equation for a combined 
lattice-harmonic potential; the single particle eigenstates 
are applied to calculate the lattice thermodynamics and 
thereby predict temperature after adiabatic transfer from 
a purely harmonic potential. In Section [TT] we discuss 
the quasimomentum distribution of a thermal gas and 
imaging of that distribution using our apparatus. In 
Section IIIII we show measurements of temperature for 
a thermal gas obtained using quasimomentum distribu- 
tions. In Section [TVl we present a numerical simulation of 
bandmapping to explain the breakdown of this method 
at high temperatures. In Section |V] we will discuss de- 
termining temperature using in-trap size and present ex- 
perimental results. 



THERMODYNAMICS OF THE COMBINED 
LATTICE-PARABOLIC POTENTIAL 
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In most optical lattice experiments to date, an ultra- 
cold gas of atoms is first created in a parabolic poten- 
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tial, followed by a slow turn-on of the periodic lattice po- 
tential. The temperature in the final combined lattice- 
parabolic potential can be predicted by measuring the 
temperature in the harmonic potential, calculating the 
corresponding entropy, and assuming that the turn-on is 
isentropic. For this procedure, the entropy in the lattice- 
parabolic potential must be calculated, which has been 
done for non- interacting particles [13, [13, [H, [IS ^30] , in 
the mean- field approximation 3l| . and for the strongly 
interacting hmit[3l, [H, [H, [35 1. In the general case, a 
quantum Monte Carlo calculation is necessary [20| . 

We employ a simple but exact non-interacting theory 
to predict temperature in the lattice-parabolic potential. 
In this section, we review the eigenstates for this poten- 
tial, which were reported in Ref. [s^. We briefly discuss 
how these states are used to calculate the entropy in the 
lattice, and thereby to predict the temperature and con- 
densate fraction after loading from a purely harmonic 
trap. 

The single-particle Hamiltonian describing this system 
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where sEb, is the lattice depth {En = [hir/dY /2m is the 
recoil energy of the atom), m is the mass of the atom, 
d is the lattice spacing, and oj is the external harmonic 
confinement frequency. Some intuition about the nature 
of the single-particle eigenstates can be gained by con- 
sidering the case when the harmonic potential is absent. 
When w = 0, the solutions to this Hamiltonian are Bloch 
wavef unctions. 
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where q is the quasimomentum of the state, / specifies 
the band-index, and u^^i{x) is a function with the same 
periodicity as the lattice. An appropriate sum over Bloch 
wavefunctions defines the Wannier function centered at 
site j [37i] . 
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Experiments with ultra-cold atoms are typically con- 
fined to the lowest band, so we drop the band index in the 
discussion that follows and let Z = 0. The Hamiltonian 
can also be conveniently expressed in second-quantized 
notation, in which the operator d\ creates a particle with 
wavefunction Wi{x). In the tight-binding approximation 
(accurate for s > 4), Eq. [1] becomes: 
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(J = ■:^ERS^/*e-^^ in the limit Vo > Er ^) is the 
nearest- neighbor tunneling matrix element, i and j index 
lattice sites, () represents a sum over nearest-neighbors, 
and fii = dlui is the number of atoms at site i. For u; — 0, 
the energy spectrum is given by the usual tight-binding 
dispersion relation. 
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where q is the quasimomentum, and \q\ < qs — hn/d 
{qs is the Brillouin-zone momentum). The dispersion 
relation is characterized by an energy gap at qB between 
the lowest energy band and the first excited band. 

The parabolic potential present in ultra-cold atom 
lattice experiments breaks the periodic symmetry of 
the lattice potential and therefore changes the single- 
particle eigenstates ^^\. Recently, an analytic solution 
was derived for the eigenstates and eigenenergies for the 
combined parabolic-lattice potential [36], which enables 
straightforward calculation of non-interacting thermody- 
namics. From these solutions, the energy of a state with 
quantum numbers n,j;,ny,nz for the Hamiltonian in Eq. 
□ is E„^ = E^^ + E,,^ + E^^ (41|] with 
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where n = 0,1,2,3,... is an integer, Q. = mtj^d^/2, 
a = 4J/r2, and a„(a) and 6„+i(a) are the Mathieu char- 
acteristic values. 

The eigenstates are ^'(f) = ^-^^ (a;)4'„^ (?/)^'„^ (z), 
where the sum ^'„^(a;) = J2j fj^'^^ji^) runs over one 
direction of the lattice. The site-dependent weights are 



i J^"" dx ce„(x,-Q;)cos(2jx), 
i /o'^da; se„+i(x, -a)sin(2ja;). 
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where ce and se are the even and odd periodic tt solu- 
tions of the Mathieu equations with parameter a. The 
probability distributions along one lattice direction for 
several energy levels are plotted in Figure [T] The low 
energy states are similar to discretized harmonic energy 
eigenstates with an oscillator frequency lo \/ m/m* , where 

TO* = 2^77 [13| is the effective mass resulting from the 
lattice dispersion relation (Eq. [6]). As the energies ex- 
ceed the bandwidth 4J of the uniform lattice potential, 
the states become localized to a few lattice sites; the lo- 
calized states are found farther from the center of the 
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1: Spatial (left column) and quasimomentum (right col- 
umn) probability distributions for a one dimensional com- 
bined lattice-parabolic potential (light grey curve in top left 
image) for different quantum numbers n and a = 352.175. 
For low n, the states are similar to discretized harmonic oscil- 
lator wavefunctions, but for higher energies the states become 
localized away from the center of the parabolic potential. 



harmonic trap as n increases. The localized states were 
observed to affect transport [l^, and were detected 
directly in Ref. [42] . The energies of these states are in- 
side the uniform lattice band-gap, reminiscent of surface 
states in solids. 

To calculate thermodynamic quantities using these 
states we directly calculated the grand canonical poten- 
tial for bosons by summing over all states in 3D j44| , 
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where /3 = l/fc^T and 3 = e"^ is the fugacity (/.t is the 
chemical potential and fc^ is Boltzmann's constant). The 
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and energy, 
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were used to calculate entropy. 
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We also used Eq. [HI [IHl and [TT] to calculate the thermo- 
dynamics for the harmonic trap by using the harmonic 
trap eigenenergies i?„^,„„,„, = h^{nx + Uy + 71^) 

To calculate the final temperature and condensate frac- 
tion in the lattice based on adiabatically loading atoms 
from the harmonic trap (at a given initial temperature 
and number of atoms), we first calculated the entropy 
and fugacity of the atoms in the harmonic trap. Then we 
varied temperature and fugacity so that the final entropy 
and number were identical in the combined parabolic- 
lattice potential. 



II. QUASIMOMENTUM DISTRIBUTION 

In this section, we determine the quasimomentum 
distribution for non-interacting particles in the joint 
parabolic-lattice potential. For particles in a periodic 
potential, the quasimomentum operator plays the role of 
the momentum operator for free particles. The Hamilto- 
nian from Eq. 3] for atoms in a lattice-parabolic potential 
written in terms of the quasimomentum and position op- 
erators is 
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To calculate the quasimomentum distribution, the 
single-particle eigenstates must be transformed into 
functions of quasimomentum |3(|: ^n^,ny,n^iq) = 
fn^{qx)fny{qy)fnAQz) (scc Appendix o), where 
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analogous equations apply for fuyiqy) and fn^iqz)- For 
the finite temperature, statistical distribution of the 
Hamiltonian in Eq. [13] we add a Bose-Einstein factor. 
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where the sum is over ah eigenstates. Fitting images 
of the quasimomentum distribution to p{q) calculated 
using Eq. [15] is infeasible given modest computational 
resources. Fortunately, by writing the Hamiltonian in 
terms of quasimomentum, we are justified in using the 
semi-classical distribution, since the rapid spatial varia- 
tions in the potential and wavefunction have disappeared 
[21I . l46j . The semi-classical approximation is appropriate 
in this system because the quasimomentum and position 
operators are conjugate [i^l 



ih. 
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In the semi-classical distribution, the operators in the 
Hamiltonian are replaced with classical variables, and 
each particle occupies the minimum uncertainty volume 
in phase space consistent with the commutation relation. 
Using this approximation, the three-dimensional, finite- 
temperature quasimomentum distribution for bosons in 
a lattice-parabolic trap is (where hk — q), 
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After integrating over the spatial degrees of freedom, 
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where Lin{u) — X]fc°=i ^''Z^" polylogarithm func- 

tion. This equation reveals that, in the thermal limit 
(3 <C 1), the quasimomentum distribution is solely deter- 
mined by the dimcnsionless parameter = ksT/J, 
which characterizes the ratio of the temperature to the 
band width. In the quantum degenerate regime (3 ~ 1), 
the fugacity and the parameter ksT/J both control the 
shape of the quasimomentum distribution. 

To image the quasimomentum distribution, we em- 
ploy "bandmapping" followed by standard time-of-flight 
(TOF) absorption imaging. Bandmapping involves turn- 
ing off the lattice potential slowly with respect to the 
bandgap, but quickly with respect to the slowly vary- 
ing harmonic trapping potential. This procedure maps a 
state of quasimomentum, q, to a state of momentum q. 
This technique has been demonstrated by several groups 
[H HE Hi], but we wih show in Sec. IV that it has hmits 
of validity not previously discussed in the literature. 

Each pixel in the two-dimensional absorption image 
is a column integral of the quasimomentum distribution 
along the imaging axis. We use an atypical imaging ge- 
ometry in our experiment: the imaging axis makes 90° 



and ±45° angles with the lattice directions, which define 
the quasimomentum axes in space after TOF. For our ex- 
periment, the coordinates ai and a2 define the imaging 
plane; the integration direction is along coordinate a^; 
and fci, /c2, and k^ are coordinates along the lattice axes 
given by, 



ki 
k2 



2^1 



2«2 + ^a3, 



2^1 + 



V2 



0:3. 



(19) 
(20) 

(21) 



The two-dimensional distribution integrated along the 
direction of the probe beam is 



/9(ai,a2) = / das /9(fci, fc2, /cs), 
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which can be expressed as a doubly infinite sum 



J 



°° -2nJf3[3-cos{b)] n I sin f-TT - n^ll i 

p(ai, 02) =AJ2 575 <! 2(^ - a)/o [cos(a)2nJ/3] + ^ 1, [cos(a)2n J/5] 4^^^^^^ ^ \ , (23) 



n=l I j=\ 



where A is a proportionality constant, a — |ai +q:2|/\/2, ^= |Q;i—a2|/\/2 and /j (a;) is the modified Bessel function 



of the first kind. 

Examples of predicted images for 3 = 0.5 and for dif- 
ferent values of ksT/J are shown in Figured A charac- 
teristic feature of these distributions is the sharp edge at 
the Brillouin zone momentum for high fc^T/J. We have 
verified that distributions calculated using a direct sum 
over eigenstates (Eq. [T5|) agree with the semi-classical 
result for the range of ksT/J used in this paper. 



III. MEASURING TEMPERATURE USING 
QUASIMOMENTUM DISTRIBUTIONS 

Images of the quasimomentum distribution obtained 
for a range of temperatures and lattice depths are fit to 
the semi-classical result (Eq. [23l) . In this section, we 
discuss how the fit is used to determine fcsT/J, which 
is compared with thermodynamic predictions for non- 
interacting particles. All data in this section and Sec. 
rvl are taken above the critical temperature for Bose- 
Einstein condensation. 

The procedure by which we prepare ultra-cold gases 
of ^'^Rb atoms confined in a lattice-parabolic potential 
has changed somewhat compared with Refs. ^^Rb 
atoms in the F = 1, mp = —1 state are evaporatively 
cooled in a two-stage process. The atom gas is first cooled 
using forced radio-frequency evaporation in a magnetic 
quadrupole trap with a 300 G/cm gradient along the 
symmetry axis. Before atom loss induced by Marojana 
transitions becomes significant, the atom gas is trans- 
ferred into a hybrid magnetic-optical trap by slowly turn- 
ing on a 7 W, 1064 nm Gaussian laser beam focused to 
a 90 /.tm waist; the magnetic field gradient is simulta- 
neously reduced to produce a force on the atoms equal 
to gravity. The second stage of evaporative cooling pro- 
gresses by increasing the quadrupole gradient and reduc- 
ing the optical power [i^ . Finally, the magnetic gradient 
is again relaxed to balance gravity and the dipole beam 
power is reduced to 0.4 W, which creates a harmonic 
trap with oscillator frequencies 56.4 ± 0.7, 29.2 ± 1.5, 
and 39.1 ± 0.2 Hz. The temperature Tho of the atom gas 
in this trap is controlled by altering the evaporative cool- 
ing sequence, both by changing the final radio-frequency 
applied during cooling in the quadrupole trap and the 
minimum optical power used during cooling in the hy- 
brid trap. 

After the dipole beam power is decreased to 0.4 W, the 
optical lattice beams (A = 812 nm) are super-imposed on 
the hybrid trap (see Ref. fd\ for more information on the 
lattice geometry). Bandmapping is implemented by turn- 
ing off the optical lattice beams in 750 ^s using a linear 
ramp of the optical lattice power. Following bandmap- 
ping, the 1064 nm beam and the quadrupole magnetic 
trap are turned off and an absorption image is taken af- 
ter 8-20 ms of time-of-flight (TOF). 

Figure [3] shows the measured temperature T in the lat- 
tice as the temperature in the harmonic trap is varied. 
Images taken after bandmapping and time-of-flight ex- 
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FIG. 2: (Color Online) Examples of predicted images from 
Eq. 26 for fesT/J = 0.5, kaT j J = 2 and kaT j J = 15 (in 
the thermal limit). Cross-sections through the distribution 
are also shown. Along cross-section A, the distribution has a 
sharp edge at the Brillouin-zone momentum at high fcsT/J. 
For low values of fcsT/J, this edge vanishes and the distribu- 
tion is similar to a momentum distribution for atoms confined 
in a harmonic trap. The direction of gravitational accelera- 
tion, relative to the imaging plane for our experiment is 
labeled by an arrow. 
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in Fig. [3] — when Tho leads to a thermodynamically pre- 
dicted ksT/J < 2.5. The theoretical prediction is calcu- 
lated using the measured Tho and N, as discussed in Sec. 
m We also include the change in the parabolic confine- 
ment induced by the Gaussian profile of the lattice laser 
beams. We assume that the overall confinement is spher- 
ically symmetric, with a harmonic oscillator frequency 
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where loq = 40Hz is the geometric mean of the oscillator 
frequencies for the hybrid trap, and w — 120 ± lO/im is 
the waist of the optical lattice beams. 

For higher Tho, the measured temperature in the lattice 
is systematically lower than the thermodynamic predic- 
tion. Effects due to interactions can be straightforwardly 
ruled out as the source of this discrepancy as the ratio of 
mean energy to estimated interaction energy per particle 
is greater than 30 for the data in Fig. [31 Also, a failure 
of adiabaticity can be ruled out as the measured temper- 
ature is less than the thermodynamic expectation. Fur- 
thermore, the temperature in the lattice measured using 
an alternate technique agrees with the thermodynamic 
prediction, as demonstrated in Sec. |Vl In the next sec- 
tion, we show that the deviation is consistent with a fail- 
ure of bandmapping for states of high quasimomentum, 
which are populated at high fc^T/ J. 



FIG. 3: Measured temperature T in the lattice obtained by 
fitting the quasimomentum distribution for s = 2 (top) and 
s = 6 (bottom) as the temperature T^o in the harmonic trap 
is varied. The vertical error bars are determined by the stan- 
dard deviation in several measurements that are averaged for 
each point and a 10% uncertainty in the Brilluoin-zone width 
after time-of-flight. The error bars in Tho represent statistical 
uncertainty. There is also a 5% systematic uncertainty in the 
TOF, which is not included in the error bars. The curve is 
the thermodynamic prediction for isentropic transfer into the 
lattice and non-interacting particles; the width of the theory 
curve reflects error in A'^ and u>. 



pansion are fit to Eq. [23] using a non- linear least squares 
solver, with A, ksT/ J, the maximum optical depth, and 
fugacity as free parameters. Data are shown for s = 2 
and 6. While Eqs. [^ and [T51 require corrections for next- 
nearest- neighbor tunneling at s = 2, the next-nearest- 
neighbor tunneling energy is 14% of J, which results only 
in a small deviation from the tight-binding dispersion re- 
lation. The fitted value of ksT/J is converted into tem- 
perature T in the lattice using a truncated Fourier series 
band calculation for J, based on a calibration of s within 
5% using Raman-Nath diffraction [49]. The number of 
atoms ranged from 3 x 10-^-1.4 x 10^ and the fitted 
fugacity from 0.35-0.74 for the data shown in Fig. [3] 

The measured temperature in the lattice agrees with 
the prediction for isentropic transfer — shown by the curve 



IV. SIMULATION OF BANDMAPPING IN ID 

To investigate the validity of bandmapping as a tech- 
nique for measuring the quasimomentum distribution, we 
numerically simulate the process using a Crank-Nicolson 
solver 50]. We calculated the time evolution of the 
single-particle eigenstates of the Hamiltonian in Eq. [1] 
in ID, with the lattice depth s given by: 

*' = {^(i-^)-™..>J>0' 

and the corresponding harmonic oscillator frequency de- 
termined by Eq. [24l The Wannier function at each lattice 
site was approximated by the ground state wavefunction 
of a harmonic oscillator with angular frequency 2Eii^/s/h 
(39| . A total of 240 lattice sites were included in the sim- 
ulation up to i = 750/xs, when the bandmapping process 
is completed (i.e., s = 0). The wavefunctions were then 
propagated for 20ms of free-evolution using the kernel, 

^' = ^^^^ 

which was convolved with the wavefunction obtained 
from the solver. The wavefunctions ^''^(x') obtained via 
this procedure are the single particle states when an im- 
age is taken. 
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FIG. 4: Numerical simulation of bandmapping for the n — 
state. The single-particle eigenstate is shown at the top left; 
the wavefunction after 20 ms of TOF, revealing the corre- 
sponding momentum distribution, is displayed at the bottom 
left. The wavefunction after bandmapping is displayed at the 
top right, with the wavefunction after 20 ms TOF shown at 
the bottom right. 



Figure |4] shows the result of this simulation for the 
n = eigenstate. The effect of bandmapping is to change 
the spatial wavefunction so that the momentum distri- 
bution after bandmapping is the quasimomentum distri- 
bution before bandmapping. Since the quasimomentum 
wavefunction is given by the Fourier series of the coef- 
ficients of the Wannier functions (see Eq. IA.4I in Ap- 
pendix p), and momentum is the Fourier transform of the 
wavefunction, bandmapping is equivalent to removing 
the modulation of the wavefunction at the lattice spac- 
ing. Standard TOF imaging, after sufficient expansion 
time, therefore reveals the quasimomentum distribution. 

To compare with the results of Sec. IIIIl we use the 
results of this calculation to create simulated images in 
the thermal limit for a range of kBT/J. We restrict the 
calculation to ID and simulate the time-evolution of the 
150 lowest energy eigenstates. The predicted image is 
calculated by summing the eigenstate probability distri- 
butions after bandmapping and time-of- flight, weighted 
by appropriate Boltzmann factors (3 ^ 1): 



pix') 



N ^ 



(27) 



Predicted images simulated using this scheme are shown 
as solid black lines in Fig. [5] for ksT/J = 2.88 and 
23; the coordinate x' after TOF is converted to quasi- 
momentum in the lattice using q = rnx'/r, where r is 
the free evolution time. The predicted images are com- 



*n(g) 



calculated in thermal limit 



P('?) = irE„e 
(dashed blue line). 

Apparent in Fig. [5] is a failure of bandmapping at high 
quasimomentum. During bandmapping, the bandgap 
at the Brilluoin zone edge shrinks, and the adiabatic 
timescale for changes in the lattice depth at high q is 
extended. Atoms with high quasimomentum therefore 
make diabatic transitions out of the Brillouin zone dur- 
ing bandmapping, leading to significant smoothing of the 
sharp edge in the image at g = ir/d for high kBT/J. Be- 
cause of this effect, images taken after bandmapping and 
time-of-flight at high ksT / J are not images of the quasi- 
momentum distribution. We find that in this simulation 
that we are not very sensitive to the choice of bandmap- 
ping time. Indeed, in the literature, bandmapping times 
of 20 Ms|3, 200 Ms[i,0,l, 1 msll and 2 ms[255 have 
been employed. 

To understand the impact of this problem on the 
method employed in Sec. IIIH simulated images for 
ksT/J =1-23 are fit to the ID semi-classical quasimo- 
mentum distribution 



p(9) = Ae-2/''-^[i-<=°''(^'3'i)l, 



(28) 



where A and [3' are free parameters. The top panel of 
Fig. [5] shows the fitted ksT' / J as ksT/J (used to create 
the distribution) is varied. The fit systematically un- 
derestimates the temperature for ksT / J > 3, implying 
that a failure of bandmapping is responsible for the dis- 
crepancy at high ksT/J evident in Fig. [3] between the 
thermodynamic prediction and the fitted temperature. 
We are unable to make direct, quantitative comparisons 
between simulations and our data because of the compu- 
tational complexity of calculating a 3D density distribu- 
tion after bandmapping, TOF, and projection onto the 
imaging plane. 



V. MEASURING TEMPERATURE USING 
IN-TRAP SIZE 

Because bandmapping fails to produce quasimomen- 
tum distributions at high temperature, we have explored 
another method for thermometry. We determine tem- 
perature using the measured size of the gas while it is 
confined in the parabolic-lattice potential. This tech- 
nique not only overcomes the limitations of bandmap- 
ping, but also has sensitivity superior to fitting the quasi- 
momentum distribution at high fc^T/J. The quasimo- 
mentum distribution is largely insensitive to tempera- 
ture once states at the edge of the Brilluoin zone are 
significantly populated. In contrast, the density profile 
always depends strongly on temperature, since the local- 
ized states — while nearly uniform in quasimomentum — 
extend to larger radii in the parabolic potential as their 
energy increases. 
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FIG. 5: (Color Online) Failure of bandmapping for high 
ksT/J. The top panels show samples of the finite temper- 
ature bandmapped distribution (solid black line), the finite- 
temperature quasimomentum distribution (dashed blue line) , 
and the semi-classical fit to the bandmapped distribution 
(dashed and dotted red line). The bottom panel shows the 
temperature T' obtained using a semi-classical fit to a ID 
simulation of imaging at temperature T after bandmapping 
and TOF. The red line has a slope of 1 and is present to guide 
the eye. 



The relation between in-trap size and temperature is 
given by the exact expression for the mean-squared width 
of the cloud in the lattice along a direction f, 



N 



■■■l*n.,„„„.(x)pe-^^™. (29) 

However, we can use the semiclassical expression, Eq. 1171 
for the width (3 < 1), 



(|f • X 



Id' 



X \r ■ x\ e 



(30) 



which gives the familiar expression (|f-xp) — ksT jmuj'^ . 
Comparing Eq. [30] to the width numerically calculated 
from Eq. [21] (Fig. [5^), the agreement is excellent as long 
as the temperature is greater than the spacing between 
the first two eigenstates given by Eq. [7] 

We cannot directly measure the in-trap density distri- 
bution using our apparatus. Therefore, we measure the 
density profile after snapping off the lattice and relatively 
short (i.e., less than 10 ms) TOF; the size of the gas in 
the trap is inferred from extrapolation. We determine 
the r.m.s. size of the gas by fitting images taken after 
expansion time i to a Gaussian profile. We then extrap- 
olate to the r.m.s. size in the trap by assuming for short 
expansion times that the gas expands with a momentum 
distribution similar to that of a gas confined in a har- 
monic potential. We therefor e fit the r.m.s. size after 
expansion time t to cr(t_) = \/ (Jq -\- At^ , leaving (Tq and 
A as free parameters [51| . The temperature of the gas 
is inferred from the fitted value of cto and the measured 
value of CO (Fig. [Bb). 

Figure[6]shows the measured temperature in the lattice 
for s = 6 as the temperature in the harmonic trap is var- 
ied. The agreement with the thermodynamic prediction 
(grey curve) is excellent, implying that the disagreement 
with this prediction evident in Fig. [3]d is due to the fail- 
ure of bandmapping at high quasimomentum. 



Conclusions 

In conclusion, we have determined the absolute tem- 
perature of atoms confined in a lattice in the thermal 
limit by employing two methods: fitting quasimomen- 
tum distributions obtained via bandmapping and mea- 
suring the in-trap size of the gas. These methods may 
be useful for verifying other methods of thermometry 
and can be applied to fermionic gases. Furthermore, de- 
termining temperature in the thermal limit will prove 
useful for future studies of unexplained transport phe- 
nomena in optical lattices at relatively high tempera- 
ture @, [13. Through numerical simulation, we have also 
demonstrated that bandmapping fails to produce accu- 
rate quasimomentum distributions at high temperature 
(or when high quasimomentum states are occupied). We 
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FIG. 6: Temperature measured using the in-trap size for 
s = 6. The curve is the theory prediction assuming adiabatic 
transfer into the lattice; the width of this curve is determined 
by uncertainty in A'^ and ui. The error bars in T reflect the 
uncertainty in the extrapolation to the in-trap size, and the 
error bars in Tho represent the uncertainty in using standard 
TOF expansion to measure the temperature of the gas before 
loading into the lattice. Inset (a) shows the exact calculation 
of the width in ID (black line) from Eq. [29] versus the semi- 
classical value s/WrJrwj^ (dashed grey line). The agreement 
is exact except at very low temperatures. Inset (b) shows a 
sample set of expansion data used to determine t he in-trap 
size. The line is a fit to the data of the form a — \/(7^^A^A^ 
where A and ctq are fit parameters. 



where j labels lattice sites, then, 



= [d^x^}{x)Y, 

'' J 


a, j d\' e*^'-^^/^ 




j '' 


■^^'^ jd^x r^ix;), 














(A.4) 



J 



Therefore, the quasimomentum-space wavefunction is 
given by the Fourier series of the Wannier function coef- 
ficients. Finally, the transformation between the quasi- 
momentum wavefunction and momentum (taking the 
Fourier transform of Eq. IA.1[) , 

f r rQB 

J J J —qB 

The momentum-space Bloch wavefunctions are 



expect this discovery to have an important impact on 
this technique, which has been applie d in many optical 
lattice experiments d 0, B [13, HE lUl • 



where wq (p) is the Fourier transform of the Wannier func- 
tion at site R = and Pj is a reciprocal lattice vector. 
Therefore, 



APPENDIX: QUASIMOMENTUM 
TRANSFORMATIONS 

The quasimomentum basis states, $^(x), are a com- 
plete set, so any wavefunction can be written in the form 
(considering only a single band), 

r r rlB 

^(x) = d\ fm^{x). (A.l) 

J J J —qs 

The function f{q) is the quasimomentum-space wave- 
function. The transformation from spatial to quasimo- 
mentum wavefunction is, 

f{q) = J d'x (j>}{x;)^ix). (A.2) 

In particular, if we have a wavefunction defined in terms 
of the Wannier functions, 

^(x) = J2^j^ji^)^ (A.3) 
j 



J J J -qs j 

(A.7) 

For \p\ < qb the only way for the integrand to be non- 
zero is if Pj — 0, since all other reciprocal lattice vectors 
lie outside the Brillouin zone, so, 

= f{p)Mp)- (A.8) 
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